Non-leptonic two-body decays of the Be meson in light-front quark model and QCD 

factorization approach 
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We study exclusive nonleptonic two-body Be (^^(s) , rjc, -l-F decays with F(pseudoscaIar or 
vector meson) factored out in the QCD factorization approach. The nonleptonic decay amplitudes 
are related to the product of meson decay constants and the form factors for semileptonic Be decays. 
As inputs in obtaining the branching ratios for a large set of nonleptonic Be decays, we use the weak 
form factors for the semileptonic Be (-D(a) , r^c, B(s)) decays in the whole kinematical region and 
the unmeasured meson decay constants obtained from our previous light-front quark model. We 
compare our results for the branching ratios with those of other theoretical studies. 



I. INTRODUCTION 

The discovery of the Be meson by the Collider Detec- 
tor at Fermilab(CDF) Collaboration [l;] in pp collisions 
at ^/s = 1.8 TeV and the subsequent measurement of 
its lifetime have provided a new window for the analysis 
of the heavy-quark dynamics and thus for an important 
test of quantum chromodynamics. Recently the CDF 
and DO Collaborations announced some new measure- 
ments of the Be meson lifetime and mass 0, 0], tb^ = 
0.463t[5 ;]^^(stat) ± 0.036(syst) ps 0, Mb, = 6275.6 ± 
2.9(stat)±2.5(syst) MeV H, and 6300±14(stat)±5(syst) 
MeV 0. The LHC is expected to produce around 
~ 5 X 10^" Be events per year d, This will provide 
more detailed information on the decay properties of the 
Be meson. Since the Be mesons carry flavor explicitly(6 
and c) and cannot annihilate into gluons, they are sta- 
ble against strong and electromagnetic annihilation pro- 
cesses. The decays of the Be meson are therefore only via 
weak interactions, which can be divided into three classes 
at the quark level: (1) the 6 — *■ g (g = c, m) transition with 
the c quark being a spectator, (2) the c — > g (g = s, c?) 
transition with the b quark being a spectator, and (3) the 
weak annihilation channels. Although the phase space of 
the c — > s, c? transitions is much smaller than the phase 
space of the b ^ c,u transitions, the Cabibbo-Kobayashi- 
Maskawa(CKM) matrix elements are greatly in favor of 
the c quark decay, hcjUcbl << |K;s|- In fact, the c-quark 
decays provide about ~ 70% of the Be decay width while 
the 5-quark decays and the weak annihilation yield about 
20% and 10%, respectively [|. This indicates that both 
6-and c-quark decay processes contribute to the Be decay 
width on a comparable footing. 

Because the b and c quarks can decay individually 
and the Be meson has a sufficiently large mass, one can 
study a great variety of decay channels. There have 
been many theoretical efforts to calculate the semilep- 
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iKiblJV'csl and \Ved \ but also the weak form factors for the 
transitions of Be to bottom and charmed mesons. The 
nonleptonic Be decays, in which only hadrons appear in 
the final state, are strongly influenced by the confining 
color forces among the quarks. While in the semilep- 
tonic transitions the long-distance QCD effects are de- 
scribed by a few hadronic form factors parametrizing the 
hadronic matrix elements of quark currents, the nonlep- 
tonic processes are complicated by the phenomenon of 
the quark rearrangement due to the exchange of soft and 
hard gluons. The theoretical description of the nonlep- 
tonic decays involves the matrix elements of the local 
four-quark operators. Although the four-quark opera- 
tors are more complicate than the current operators in- 
volved in the semileptonic decays, the nonleptonic decays 
of the heavy mesons are useful for exploring the most 
interesting aspect of the QCD, i.e. its nonpcrturbative 
long-range character. 

In our recent paper [isj , we analyzed the semilep- 
tonic Be decays such as Be {D,ric, B, Bs)£i^e and 
rji, — > Beivi{£ = e,fi,T) using our light-front quark 
model(LFQM') based on the QCD-motivated effective LF 
Hamiltonian 0, El, [H, 113, IH, |4|. The weak form 
factors f±{q^) for the semileptonic decays between two 
pseudoscalar mesons are obtained in the = frame 



-<Yi < 0) and then analytically continued to the 



MIIS 37, 38, 39, 40, ,41, 43] decays of the Be meson. 
The semileptonic Be decays provide a good opportunity 
to measure not only the CKM elements such as lUchl, 



timelike region by changing to —q in the form factor. 
The covariance (i.e., frame independence) of our model 
has been checked by performing the LF calculation in the 
q'^ — Q frame in parallel with the manifestly covariant 
calculation using the exactly solvable covariant fermion 
field theory model in (3 -I- 1) dimensions. We also found 
the zero-mode contribution to the form factor f-{q^) and 
identified the zero-mode operator that is convoluted with 
the initial and final state LF wave functions. 

In this paper, we extend our previous LFQM analysis 
of the semileptonic Be decays |43| to the nonleptonic two- 
body decays of Be mesons such as Be (^(s) , Vc, ^(s))^ 
and Be (D, 77c, i3(s) )U(here P and V denote pseu- 
doscalar and vector mesons, respectively). The QCD fac- 
torization approach is widely used since it works reason- 
ably well in heavy-quark physics [13, [Ml, [H, [H, [13, . 
The factorization approximates the complicated non- 
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leptonic decay amplitude into the product of the meson 
decay constant and the form factor. A justification of this 
assumption is usually based on the idea of color trans- 
parency [5^. We shall use the form factors for semilep- 
tonic Be — > (-D(s), 77c, decays as well as the meson 
decay constants obtained in our LFQM [H, El] as input 
parameters for the nonleptonic Be decays. As done by 
many others [5, 6, 7, 8, 9, 0, d E [ll, H, we con- 
sider only the contribution of current-current operators 
at the tree level and calculate the decay widths for various 
nonleptonic Be decays. As far as the decay width is con- 
cerned, the contribution from the tree diagram is much 
larger than that from the penguin diagram. The penguin 
contribution may be important in evaluating the CP vi- 
olation and looking for new physics beyond the standard 
model, which we do not consider in this work. 

The paper is organized as follows. In Sec. II, we discuss 
the weak Hamiltonian responsible for the nonleptonic 
two-body decays of the Be meson. In Sec. Ill, we present 
the input parameters such as the weak decay constants 
and the form factors obtained in our LFQM |43l. [48| based 
on the QCD-motivated effective Hamiltonian |44l. [isj . 
The mixing angles between rj and 77' mesons are also an- 
alyzed, both in octet-singlet and quark-flavor bases, to 
extract the decay constants relevant to 77 and 77' mesons. 
Section IV is devoted to the numerical results. A sum- 
mary and conclusions follow in Sec.V. 



II. NONLEPTONIC TWO-BODY DECAYS OF 
THE Be MESON 

The nonleptonic weak decays are described in the stan- 
dard model by a single W boson exchange diagram at tree 
level. In the standard model, the nonleptonic Be decays 
are described by the effective Hamiltonian, which was 
obtained by integrating out the heavy W boson and top 
quark. For the case of & — > c, m and c ^ s,d transitions 
at the quark level, neglecting QCD penguin operators, 
one gets the following effective weak Hamiltonian: 



TABLE I: Values for CKM matrix elements used in this work. 
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6— >c('u) 
eff 



= K,[ci(A7)Of + C2(m)0§''] 



+Vub[ci{fi)Of + C2{fi)Of] + h.C. I (1) where 



and 



K 



cff 



s(d) _ 



^(^Ve,Mt,)Of+C2{t,)Of] 
+ Ves[ci{fi)Or + C2(m)021 + h.C. L (2) 



where Gf is the Fermi coupling constant and Vq-^q^ 
the corresponding CKM matrix elements. We use the 
central values of the CKM matrix; elements quoted by 
the Particle Data Group(PDG) [57] that we summarize 
in Table [D The effective weak Hamiltonian consists of 











Veb 


Kb 


0.974 


0.2255 


-0.230 


1.04 


0.0412 


0.00393 



products of local four-quark operators ©1^2 renormalized 
at the scale fj,, and scale-dependent Wilson coefficients 
Ci,2(m)i which incorporate the short-distance effects aris- 
ing from the renormalization of Tieff from fi = to 
/i = 0{j7ib). The local four-quark operators Oi and O2 
due to b and c decays are given by 



qb 



qb 



01 



QCq 



{qb)v-A[{d'u)v-A + {s'c)v-A]{q = c, u), 
{qu)v~A{d'b)v-A + {qc)v-A{s'b)v-A{q = c,u), 
{cq)v-A{d'u)v-A{q = s, d), 

{cu)v-A{d'q)v^A{q = s, d), (3) 



where {qq)v-A = 97/^(1 "Ts)? f^nd the rotated antiquark 
fields are given by 



d' = Vudd + VusS, 

S' = Vedd+VesS. 



(4) 



Without strong-interaction effects, one would have ci — 
1 and C2 = 0. However, this simple result is modified 
by gluon exchange: i.e., the original weak vertices get 
renormalized and the new types of interactions(such as 
the operators O2) are induced [5l|. In these decays, the 
final hadrons are produced in the form of pointlike color- 
singlet objects with a large relative momentum. Thus, 
the hadronization of the decay products occurs after they 
separate far away from each other. This provides the 
possibility to avoid the final state interaction. A more 
general treatment of factorization was presented in [ssl . 

For the operators Oi = {qiq2)v-A{q'iq2)v-A and 
C2 = {qiq2)v~A{q' iq2)v-A, using the Fierz transforma- 
tion under which V — A currents remain V — A currents, 
one gets the following equivalent forms 



ciOi + C2O2 = aid + C2d = a2d + cid. 



(5) 



aiifi) ^ ci{fi) + ^C2{n), a2(M) = C2(/i) + -^ci(/i), (6) 

and Ne is the number of colors. The terms 
Oi = {qiT^q'^)v^A{q'iT''q2)v-A and 62 = 
{qiT''q2)v~A{q\T''q2)v^A with SU(3) color gener- 
ators T° are the nonfactorizable color-octet current 
operators, which are neglected in the factorization 
assumption. A detailed analysis of 1/A^c corrections to 
the coefficients ai , 02 as well as the role of color-octet 
current operators in B decays can be found in [5l| . 

In the factorization approach to nonleptonic meson de- 
cays, one can distinguish three classes of decays for which 
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the amplitudes have the following general structure [50|: 



Gf 

(class I) : -j=VcKMai{fJ,){Oi)F, (7) 

(class II) : ^VcKMa2{f^){02)F, (8) 
G 

(class III) : -^VcKMlaiin) + xa2ifi)]{Oi)F, (9) 

where {Oi)F represents the hadronic matrix element 
given as the products of matrix elements of quark cur- 
rents and cc is a nonperturbative factor equal to unity in 
the flavor symmetry limit [Flj . The first (second) class is 
caused by a color-favored (color-suppressed) tree diagram 
and contains those decays in which only a charged (neu- 
tral) meson can be generated directly from a color-singlet 
current. The first and second class decay amplitudes are 
proportional to ai and a2, respectively. The third class 
of transitions consists of those decays in which both ai 
and 02 amplitudes interfere. 

In this paper, we consider the following type of non- 
leptonic Be Fi + F2 where Fi is the pseudoscalar 
meson(77c, ^(s)7 and F2 the meson(vcctor or pseu- 

doscalar) being factored out. For instance, the factorized 



matrix element of 
is defined as 



FiF+ with Fi = f]c and F2 



XiB+FuF+) ^ {Fi\icb)v^A\B,){F2\idu)v-A\0). (10) 

The matrix elements of the semileptonic Be — > Fi de- 
cays can be parametrized by two Lorentz-invariant form 
factors: 

(Fi(P2)|(9V)y-A|i?c(Pi)) = .f+(<z')P^ + ./-(g')g^ (11) 

where P = Pi + P2 and q = Pi — P2- The two form 
factors also satisfy the following relation: 



Mi 



(12) 



The decay constants fp and fv of pseudoscalar (P) and 
vector(t^) mesons are defined by 

{P{pmq')v-A\0) = -ifpp'^, 

{V{pM{m')v-A\Q) - fvMve^ih), (13) 

where e{h) is the polarization vector of the vector meson. 
In the above definitions for the decay constants, the ex- 
perimental values of pion and rho meson decay constants 
are /tt ~ 131 MeV from n ^ and fp w 220 MeV from 
p — > e+e~. 

Using Eqs. pT|) and we obtain the following ex- 
pressions for the factorized matrix elements 



(14) 



when F2 is a pseudoscalar meson, and 

XiBtF.y^) ^ 2fyMy{e ■ J/f =^"^^^4), (15) 




(a) 




(b) 



FIG. 1: Quark diagrams for the nonleptonic Be BK de- 
cays: (a) color-favored(class I) 5+ — > B^'K^ and (b) color- 
suppressed(class II) B^ B^ K'^ decays. 



when F2 is a vector meson. For the latter case, only 
longitudinally polarized vector mesons are produced in 
the rest frame of the decaying Be meson, i.e.. 



(16) 



where 



[Ml^ - {Ml + M2)2][M2^ - (Ml - M2)2] 

2M^^ ' 

(17) 

is the center of mass momentum of the final state 
meson(Fi with Mi or F2 with A/2)- 

The decay rate for Be Fi + F2 in the rest frame of 
the Be meson is given by 



r{Be ^ F1F2) = 



Pc 



\{FlF2\n,ff\B^e 



+ \|2 



(18) 



In Fig.[l] we show the example of quark diagrams for the 
nonleptonic Be BK decays: (a) color- favored (class I) 

B^K^ and (b) color-suppressed(class II) 5+ 
B'^K^ decays. We also show in Fig. [2] the example 
of quark diagrams for the class III transitions such as 
Be — > D^D^ . In the factorization approximation, these 
nonleptonic decay amplitudes can be expressed as the 
product of one-particle matrix elements. 



A. Class I decay modes 



(1) For the h — * (u, c){qiq2) process, 
{DHl+\n.^\Bt) ^ ^K.V^,%,aiX(^^^"-*^"), (19) 



and 



(77cA/+|Heff|B+) - %KbV;* aiX(^^"-^^+), (20) 



x/2 
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C. Class III decay modes 




(a) 




(b) 



FIG. 2: Quark diagrams for the class III nonleptonic Be — > 
D'^ decay, which consist of both color-favored (a) and 
color-suppressed (b) decays. 



(2) For the c — * [d, s)(qig2) process, 



(B°M+|Heff|S+) = ^K.V,%,aiX(^^^°'^^"), (21) 



and 



(i?."M+|Heff|S+> - ^K.V^,%,aiX(^?^".^'^"), (22) 

where M = ir, K, p, K* and Vg-^q^ = V^d or V^s depending 
on whether M — (tt, p) or (if, if*). 



B. Class II decay modes 



(1) For the {d,s){qiq2) process. 



{D+Am.^\Bt) = ^y,,,^;,a2X(^?^^*^"), (23) 



and 



B+M'\n.^\Bt) = i^K,,K;.«2X(^^^"'^^"), (25) 



(2) For the c — > u{qiq2) process 
G 
72 
and 

-l-K.K*.^'''^^''^'"^"^], (26) 

where M ~ (tt, 77^'^ p, w, £>(*)) for c id,s) induced de- 
cays and M = (tt, p, oj, K, K, K*, K*) for c — > u induced 
decays, respectively. As in the case of color-favored class 
I decay modes, the factorized matrix elements for color- 
suppressed class II decay modes can be obtained from 
Eqs. p4| and (jlSp except that the decay constants for the 
neutral 7r°, p'^, and uj mesons are replaced by /p(v)/V2. 



For the class III _B+ D^Af^q — d,s) transitions, 
the decay amplitude is given by 



V2 



-a2X(^="^?^^'^")], (27) 



where Vp = Kh for M = rjc and Vub for M = D. The 
expressions for the factorized matrix elements can be ob- 
tained from Eqs. ([11]) and p^ . 



III. INPUT PARAMETERS 

In this section we shall briefly discuss and summarize 
all of the input parameters, such as the model parame- 
ters, decay constants, and form factors for semileptonic 
Be (-D(s), 77c, decays, which are relevant to the 

present work. 



A. Brief review of LFQM 



The key idea in our LFQM [4J,|45l for the ground state 
mesons is to treat the radial wave function as a trial func- 
tion for the variational principle to the QCD-motivated 
effective Hamiltonian saturating the Fock state expan- 
sion by the constituent quark and antiquark. The QCD- 
motivated effective Hamiltonian for a description of the 
ground state meson mass spectra is given by 



Hqq — Hq 



where 



+ Vqq = ^Jml + k^ + ^m} + k^ + Vqq, (28) 



4a. 2 
AT 6 mq-rriq 



(29) 



In this work, we use two interaction potentials: (1) the 
Coulomb plus linear confining (i.e. n = 1) potential 
and (2) the Coulomb plus harmonic oscillator(HO)(i.e. 
n = 2) potential, together with the hyperfine interaction 
(Sg • Sq) ~ I/4(— 3/4) for the vector (pseudoscalar) me- 
son, which enables us to analyze the meson mass spectra 
and various wave- function-related observables, such as 
decay constants, electromagnetic form factors of mesons 
in a spacelike region, and the weak form factors for the 
exclusive semileptonic and rare decays of pseudoscalar 
mesons in the timelike region [H, [H, [^i, [43, IH, . 

The momentum-space light-front wave function of the 
ground state pseudoscalar and vector mesons is given 
by -^{xi.'ki^^Xi) = nx^x^{xi,'ki±)(l){xi,'ki±), where 
(/)(xi,kij^) is the radial wave function and TZx^x^ is the 
covariant spin-orbit wave function. The model wave 
function is represented by the Lorentz-invariant vari- 
ables, Xi = pf /P^, \ii± = pi± — XiP± and Xi, where 
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p/. = (P+,P-, Pj^) = (P0 + P3,(M2 + P2_)/P+,P_L) is 

the momentum of the meson M, and and are the 
momenta and the hehcities of constituent quarks, respec- 
tively. 

The covariant forms of the spin-orbit wave functions 
for pseudoscalar and vector mesons are given by 



-^iAl(pl)75^'A2(P2) 



Ai Ao 



-"Ai(Pl) 



Mo+mi+m2 



V2M0 



(30) 



where Mq = y/M^ - (mi - 7712)^ M§ = + 
m1)/xi is the boost invariant meson mass square ob- 
tained from the free energies of the constituents in 
mesons, and e^(Jz) is the polarization vector of the vec- 
tor meson [GOj . For the radial wave function (p, we use 
the same Gaussian wave function for both pseudoscalar 
and vector mesons: 



47r3/4 fdk 



^3/2 



:gexp(-fcV2/32), (31) 



where /3 is the variational parameter and dkz/dx is the 
Jacobian of the variable transformation {x,kj^} k = 
(ki,A:,). 

We apply our variational principle to the QCD- 
motivated effective Hamiltonian first to evaluate the ex- 
pectation value of the central Hamiltonian Hq + Vq, i.e., 
{4>\{Hq + Vo)|(/'), with a trial function 0(xi, kij_) that de- 
pends on the variational parameter [3. Once the model 
parameters are fixed by minimizing the expectation value 
((/>|(-ffo + K))|0), the mass eigenvalue of each meson is 
obtained as M^q = {4>\{Hq + Vqq)\(j)). Minimizing en- 
ergies with respect to fi and searching for a fit to the 
observed ground state meson spectra, our central poten- 
tial Vb obtained from our optimized potential parameters 
(a = -0.72 GeV, h = 0.18 GeV^, and a, = 0.31) AM, for 
the Coulomb plus linear potential was found to be quite 
comparable with the quark potential model suggested by 
Scora and Isgur [6l|, where they obtained a = —0.81 



GeV, h = 0.18 GeV^, and as = 0.3 - 0.6 for the Coulomb 
plus linear confining potential. A more detailed pro- 
cedure for determining the model parameters of light- 
and heavy-quark sectors can be found in our previous 
works [44l4a|. 

Our model parameters (wg, Pqq) obtained from the lin- 
ear and HO potential models are summarized in Table [TTl 
The predictions of the ground state meson mass spectra 
including bottom-charmed mesons can be found in our 
recent work, Ref. lisl . 



B. Form factors for semileptonic Be P decays 

For the nonleptonic two-body Be decays, we use the 
dependent form factors f+{q^) and fo{q^) for the Be 
{D, Ds,r]c, B, Bs) decays as input parameters. 

Within the framework of LF quantization, while the 
form factor f+{q^) can be obtained only from the valence 
contribution in the q^ — Q frame with the compo- 
nent of the currents without encountering the zero-mode 
contribution [g^l, the form factor f-{q^) [or equivalently 
/o((7^)] receives the higher Fock state contribution (i.e., 
the zero-mode in the g"*" = frame or the nonvalence 
contribution in the g"*" > frame). In order to calculate 
f-{q^), we developed in [4^,|43| an effective treatment of 
handling the higher Fock state (or nonvalence) contribu- 
tion to /-(g^) in the purely longitudinal g+ > frame 
(i.e., g^ — q^q~ > 0) based on the Bethe-Salpeter(BS) 
formalism. In our recent LFQM analysis [43| of the 
semileptonic Be {D,ric, B, Bs)£i^e decays, we utilized 
our effective method [46] to express the zero-mode contri- 
bution as a convolution of zero-mode operator with the 
initial and final state LF wave functions. In this way, 
we obtained the form factor /-(g^) in the g+ = frame 
using the perpendicular components of the currents and 
discussed the LF covariance of /-(g^) in the valence re- 
gion by analyzing the covariant BS model and the LF 
covariant analysis described by Jaus [63j . 

The LF covariant form factors /+(g^) and /-(g^) for 
-Sc(gig) P{Q2q) transitions obtained from the g"*" = 
frame are given by (see (43j for more detailed derivations) 



Ml') 



dx 



d^k_L (l)i{x,'ki_) (j)2{x,'k'^_) 
16^3 ^^2 + ^Ai+ K 



(A^2 + ki-k'_L), 



J 167r3 ^Ai + ki ^Ai + ^ 1 « V 2 q) I 



+2 



g.P 



,(k±-q±) 



, (k± • qj 



^^i^[M| - (1 - x){q^ +q-P) + 2xMl 



^(1 - 2x)Ml - 2{mi - mq){mi + ma)] \, 



(32) 
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TABLE II: The constituent quark mass[GeV] and the Gaussian parameters /3[GeV] for the hnear and HO potentials obtained 
by the variational principle, q = u and d. 



Model 


niq 


nis 


rric 


nib 






(5ss 






Pec 




Psb 




Pbb 


Linear 


0.22 


0.45 


l.S 


5.2 


0.3659 


0.3886 


0.4128 


0.4679 


0.5016 


0.6509 


0.5266 


0.5712 


0.8068 


1.1452 


HO 


0.25 


0.48 


1.8 


5.2 


0.3194 


0.3419 


0.3681 


0.4216 


0.4686 


0.6998 


0.4960 


0.5740 


1.0350 


1.8025 



where k'^ = + (1 — .T)q_L, -Ai — {1 — x)mi + xniq 
{i ^ 1,2), and q • P = - M| with Mi and Ma 
being the physical masses of the initial and final state 
mesons, respectively. We should note that the LF co- 
variant form factor f^{q^) in Eq. ([32|) is the sum of the 
valence contribution fl^^iq^) and the zero- mode contri- 
bution /ZM-(q2)^ 

For the analysis of the nonleptonic Be — > DgF de- 
cays where F is the vector or pseudoscalar meson be- 
ing factored out, we show in Fig. [3] the g^-dependence 
of the weak form factors f+{q^) (solid hne) and fo{q^) 
(dashed line) for the Be Ds transition obtained from 
the linear (upper panel) and HO (lower panel) potential 
parameters. The circles represent the valence contribu- 
tion /™'(g^) to fo{q^)- That is, the difference between 
fo{q^) and /™'(9^) represents the zero- mode contribu- 
tion to fo{q^). We obtain /+(0) = /o(0) = 0.120 [0.126] 
a.t q^ — for the linear [HO] potential model. The form 
factors at the zero-recoil point (i.e. q^ = <Zniax) ob- 
tained as /+(gn,ax) = 0-992 [0.868] and /o(g^ax) = 0-475 
[0.493] for the linear [HO] potential model. On the other 
hand, the valence contribution to fo{q^) at the zero- recoil 
point is obtained as /^H^max) = 0-442 [0.443] for the lin- 
ear [HO] potential model. 

In Table IIIIl we show the decay form factors /+ (0) = 
/o(0) at = for the semileptonic Be — > [D, rjc, B, Bg) 
decays obtained from [i^ and the rare Be Dg decay 
obtained in the present work (i.e. Fig. [3]) and compare 
them to other theoretical model predictions. 



C. Weak decay constants of 77 and rj' 

In this work, we shall also consider the nonleptonic 
decays of Be mesons to isoscalar states such as ui and 
{r],rj'). Isoscalar states with the same J^*^ will mix, but 
mixing between the two light-quark isoscalar mesons, and 
the much heavier charmonium or bottomonium states 
is generally assumed to be negligible. Since the vector 
mixing angle is known to be very close to ideal mixing, 
we assume ideal mixing between w and </) mesons, i.e., 
uj = {uu+dd)/^/2 and (j) = ss. However, the octet-singlet 
mixing angle of 77 and 77' is known to be in the range of 
— 10° to —23°. The physical 77 and 77' are the mixtures of 
the flavor SU{3) octet 778 and singlet 770 states: 

(;,)=.(.,(-), (33, 



"> DJLine 


ar potential) 






ft- 

- — fo" 
" fo"' 

























O I ^ ^ ^ ^ \ ^ ^ ^ ^ \ ^ ^ ^ ^ L 

O 5 10 15 

q^fGeV^] 



> D iHO ptnetmat) 




A 

- — fo" 




" fo" 













O I , , , , \ , , , , \ , , , , L 

O 5 10 15 



FIG. 3: The weak form factors /+(g^) (solid line) and fo{q^) 
(dashed line) for the Be Da transition obtained from the 
linear (upper panel) and HO (lower panel) potential parame- 
ters. The circles represent the valence contributions /o^'(g^) 
to /o(g2). 

where 

---,„, / cos^ — sinS \ 

UiO) ^ { ■ a a ] , 34) 

^ ' \^ smt' cost' J ^ ' 

and 778 = {uu + dd — 2ss) / and rjo = {uu+dd+ ss)/\/3. 
Analogously, in terms of the quark-flavor basis rjq = {uu+ 
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TABLE III: Form factors /+(0) = /o(0) at = for Be ^ {D(^^),r]c, B, Bs) transitions. 





Linear [HO] 


f8, 9] 


[19] 


L22] 


fl4] 


[26] 


[23] 


(16] 




^-"(0) 0.086[0.079] 


0.14 


0.69 


0.1446 




0.089 


0.16 


0.08 ±0.02 




"""=(0) 0.120(0.126] 












0.28 


0.15 ±0.02 




~*''=(0) 0.482(0.546] 


0.47 


0.76 


0.5359 


0.49 


0.622 


0.61 


0.58 




^""(0) 0.467(0.426] 


0.39 


0.58 


0.4504 


0.39 


0.362 


0.63 


0.41 ±0.04 




"""=(0) 0.573(0.571] 


0.50 


0.61 


0.5917 


0.58 


0.564 


0.73 


0.55 ±0.03 



dd)/y/2 and ry^ = ss, one obtains [64 



nq 



(35) 



The two schemes are equivalent to each other by = 
9 + arctan-\/2 when SU/(3) symmetry is perfect. How- 
ever, when one takes into account the SU/(3) breaking 
effect, this relationship is not maintained but given by the 
following Fock decompositions of the octet-singlet basis 
states fei]: 



1^8) = 



+ \uu + dd~ 2ss) 

3 

V2(*g - ^'s) \uu + dd + ss) 

^ 3 ' 

\/2(^'q - -^s) \uu + dd- 2ss) 

3 7! 

+ \uu + dd + ss) 

^ 3 v! ' 



(36) 



where denote LF wave functions of the corresponding 
parton states. Only in the SU/(3) symmetry limit, i.e., 
= "ifs, would one find pure octet and singlet states in 
Eq. ([55]) . Although it was frequently assumed that the 
decay constants follow the same pattern of state mixing, 
the mixing properties of the decay constants will gener- 
ally be different from the mixing properties of the me- 
son state since the decay constants only probe the short- 
distance properties of the valence Fock states while the 
state mixing refers to the mixing of the overall wave func- 
tion [el]. 

Using the decay constants of Eq. (fT5)) defined in the 
quark- flavor basis, the two basic decay constants fq and 
fs arising from rjq and rjs are obtained as 



2V6 / dx 



.1 167r3 



q{s)i 



(37) 



and simply follow the pattern of state mixing due to the 
Okubo-Zweig-Iizuka(OZI) rule ^H, i.e. 



(38) 



The OZI rule implies that the difference between the two 
mixing angles 4>q and 0^ vanishes (i.e., 4>q = (jig = (f) in 
Eq. dSSl) to leading order in the 1 /N^ expansion. On the 



other hand, the deca y co nstants in the octet-singlet basis 
are parametrized as |63 . [66| 



J ri J n 
J Tj' •> iY 



COS ^8 

sin 08 



- sm fo 
cos 9o 



fs 
fo 



(39) 



where 6g, and 6q turn out to differ considerably and be- 
come equal only in the SU/(3) symmetry limit. 

By using the correlation between the quark-flavor mix- 
ing scheme and the octet-singlet scheme [U, [6^ , one ob- 
tains 



/8 



arctan(V2/s//g), 



/o = 0o = 0-arctan(V2/,//,). (40) 

In our previous work [i^ ]. we obtained the rj ~ rj' mix- 
ing angle 9 ~ —19° for both linear and HO potential 
models by fitting the physical masses of 77 and 77'. This 
corresponds to the mixing angle ^ = 35.7° in the quark- 
flavor basis. We applied this mixing angle to predict the 
decay widths for 77(77') — > 77 using the axial anomaly 
plus partial conservation of the axial vector current 
(PCAC) relations ^ and obtained /g//^ = 1.32 (1.25) 
and fol f-n = 1-16 (1-13) for the linear (HO) potential 
model ,44j . From the decay constants of octet and singlet 
mesons together with Eq. ([^0]) . we now obtain the four 
parameters (/,, fs,98,0o) as follows: fq/f„ = 0.97 (1.00), 
= 1.46 (1.36), Os = -29.2° (-26.8°), and = 
— 7.3° (—10.6°) for the linear (HO) potential model, re- 
spectively. Given this background, we finally obtain the 
decay constants related to the t] and 77' mesons as follows 



/« = 103.2 (106.4) MeV, 
f^y = 74.2 (76.4) MeV, 



= -116.6 (-104.0) MeV, 
= 155.3 (144.7) MeV. 

(41) 



Our results for the mixing parameters of 77 and rj' 
are consistent with those obtained from Feldmann et 
al. [H, namely, /g//. - 1.26, /o//. = 1.17, = 
(1.07 ± 0.02), = (1.34 ± 0.06), ^ = (39.3° ± 1.0°), 

Og, = -21.2°, 6*0 = -9.2°, /9 = (108.5 ± 2.0) MeV, 
= -(111.2 ± 5.0) MeV, = (88.8 ± 1.7) MeV, and 
f^, = (135.8 ±6.1) MeV. 

For the measured values of meson decay constants, we 
use the central values extracted from the experimental 
measurements [stI l68l. l69i]. However, for the unmeasured 
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decay constants, we use the average values obtained from 
our linear and HO model predictions in [i^ in addition to 
the present work. The values for the decay constants used 
in this work are compiled in Table IIVI Since the decay 
constant /,,^ extracted from CLEO Collaboration [t^I has 
a large error bar, i.e., /^^^° — 335 ±75 MeV, we instead 
take the average value fj^^ = 340 MeV of our LFQM 
predictions gl], i.e., /^^ = 326 MeV and /^^ = 354 
MeV. 



IV. NUMERICAL RESULTS 

In our numerical calculations of exclusive Be decays, 
we use two sets of model parameters (to, (3) for the lin- 
ear and HO confining potentials given in Table [TT] to 
compute the weak form factors for semileptonic Be — > 
{D(^s),'']c, B(^s^) decays and the unmeasured decay con- 
stants as given in Tables IIIII and llVi respectively. Us- 
ing them together with the CKM matrix elements given 
by Table [T] we finally predict the branching ratios which 
are given in Tables |V] and IVII Although we show the 
form factors in Table IIIII only at a maximum recoil 
point = 0, we use the form factor at = Mp 
obtained from (43j for the corresponding nonleptonic 
Be -> {D(s),-nc,B(s))MF decays. 

In Table |V] we show the nonleptonic decay widths of 
the Be meson for a general value of the Wilson coefficients 
ai and 02, whereas in Table IVTl we give the correspond- 
ing branching ratios (in %) at the fixed choice of Wil- 
son coefBcients [1]: {a\) = 1.20 (1.14) and {a\) = 
—0.317 (—0.20) relevant for the nonleptonic decays of the 
c (5) quark. For the lifetime of the Be, we take the central 
value T{Be) = 0.46 ps (i.e. Ftot = 1.43 x lO^^^ Qgy) pre- 
sented in PDG [57]. Our branching ratios for both b and 
c induced decays listed in TableJ\Tlare generally close to 
the other quark model results 0, H, H, 13, but differ 
substantially from the ones obtained by Refs. [5l.[l0l[40|. 

The relative size of the branching ratios for various de- 
cay modes may be estimated from power counting of the 
Wilson coefficients and the CKM factors with respect 
to the small parameter of the Cabibbo angle A = sin 9c in 
the Wolfenstein parametrization [7l|, e.g. the CKM ma- 
trix elements can be expanded in terms of A as Vud ^ 1, 

Vus ~ A, Vub ~ A3, Ved - -A, Ves - 1, and Veb - A2. 

From Tables |V] and IVII we make the following observa- 
tions: 

(1) The class I decay modes determined by the Wilson 
coefficient ai have comparatively large branching ratios. 
The CKM favored c decays such as B+ B° (7r+,p+) 
decays with the CKM factor VcsV*^^ ^ A'^ have branching 
ratios of the order of 10~^ (e.g. 2% ^ 4% in our model 
predictions), which are the most promising class I decay 
modes shown in Tables IVl and IVII The CKM-suppressed 
c decays such as B^ — > B^^K^ with VesV*^ ~ A^ and 
B+ B° (7r+,p+) with VedV*^ - A\ as well as the 
CKM-suppressed b decays such as B^ rje (7r+,p+) 
with VcbV*^ ^ A^, have branching ratios of the order 



of 10~3 and should still be accessible at high luminosity 
hadron colliders. However, the branching ratios of the 
b ~^ u induced decay modes are too small 0(10^*^ — 10^^) 
to be measured experimentally. 

(2) The branching ratios for the class II decay modes 
determined by 02 are relatively smaller than those for the 
class I decay modes. However, the i?+ — > B^ {K^,K*^) 
decays with VesV*^ ~ A° have branching ratios of the or- 
der of 10~'^ and these modes should be accessible exper- 
imentally. Of interest is the abnormally small branching 
ratio of B t — > B'^-q' compared to that of B'^ — > B^rj. As 
stated in [40|, the reason for such a small branching ra- 
tio is not only because the available physical phase space 
is too small but also because there are large destructive 
interferences between rj'^ and rj'^ due to the serious can- 
cellation between the CKM factors VcdV*^ and VesV*^. 

(3) The class HI decay modes involve the Pauli inter- 
ference. Taking into account the negative value of 02 
with respect to ai, one can see that the class HI de- 
cay modes shown in Table |V] should be suppressed in 
comparison with the cases in which the interference is 
switched off. In order to test the effects of the interfer- 
ence, one may put the widths in the form of F = Fo-t-AF, 
where Fq = xiof -|-a;2a2 and AF = 2:0102, and then com- 
pute the AF/Fo (in %) as done in 0]. Our absolute 
values of AF/Fq obtained from the linear (HO) model 
are 34.0 (34.0)% for B+ D+D°, 37.3 (42.4)% for 
B+ D+D°, 18.4 (15.1)% for B+ r]cD+ , and 20.2 
(18.9)% for 5+ VcD^, respectively. This indicates 
that the interference is the most significantly involved in 
the _B+ BfD^ decay compared to others. In partic- 
ular, the B+ D+D^ and B+_-^ D+Dq decay modes 
have been proposed in f2!^, IH, [H, |4l| for the extraction 
of the CKM angle 7 through amplitude relations. 



V. SUMMARY 

In this work, we have studied the exclusive nonlep- 
tonic Be — > (-D(s) , 77c, B(^g))M decays where the final state 
M mesons are factored out in the QCD factorization 
approach. The inputs used to obtain their branching 
ratios were the weak form factors for the semileptonic 
Be — *■ (l?(s), 77c, B(s)) decays in the whole kinematical 
region and the unmeasured weak decay constants ob- 
tained from our previous LFQM analysis [H, [H, IH, \^ . 
For the measured values of decay constants, we use the 
central values extracted from the experimental measure- 
ments [Hii,!!^. 

Our predictions for the branching ratios are summa- 
rized in Tables |V] and IVII and compared with other the- 
oretical results. Overall, the class II decay modes have 
more discrepancies among the theoretical models than 
the class I and HI decay modes do. The upcoming 
experimental measurements of the corresponding decay 
rates can examine various theoretical approaches. The 
most promising measurable decay modes appear to be 
the CKM favored c decays such as i?^ — > B^ {■k'^ , p'^) 
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TABLE IV: Meson decay constants(in unit of MeV) used in this work. 

/vr /k fp fx' f^, f^i fp fp' /t7e fp^ 

131 [57] 159.8 [57] 220 [57] 195 [57] 217 [57] 104.8 -110.3 75.3 150.0 222.6 [68] 241 [48] 340 [48] 259.5 [69] 



TABLE V: Exclusive nonleptonic decay widths F (in 10 GeV) of the Be meson for the general values of the Wilson coefficients 
a\ and 02. 



Class Mode 


Lin [HO] 


[8, 9] 


[14] 


[26] 


BJ D Tv^ 


4.7[4.7](xl0 )al 


- 


- 


- 


B+ ^ D"p+ 


1.4[1.2](xl0"^)a? 


- 


- 


- 


BJ D 


3.9[3.4J(xl0 )af 








Be D K ^ 


7.5[6.4](xl0-^)af 








n+ + 

B2 r?c7r+ 


Q.997[1.280Jai 


0.93ai 


1.02af 


1.47a? 




2.827[3.563]a? 


2.3af 


2.6O0? 


3.35a? 


Bt rieK+ 


0.081[0.103]a? 


0.073a? 


0.082a? 


0.15a? 


T T~> -1- * -1- 

I BJ ^ rjcK ^ 


0.147|0.184Jai 


0.12ai 


0.15a? 


0.24a? 


BJ B 7v^ 


1.557(1. 296Jai 


l.Oaf 


1.10a? 


1.51a? 


Be B 


1.936[1.505]a? 


1.3af 


1.41a? 


1.93o? 


Bc B 


0.126[0.104]af 


0.09ai 


0.098a? 


~ 


Be — ^ B 


0.042[0.032]af 


0.04ai 


0.038a? 




Ti4- TiO -\- 

B+ B^7r+ 


36.97(36. 71]a? 


nr 2 

25ai 


34.7a? 


34.78a? 


BJ B"p^ 


25. 43(23. 22Jaf 


1 A 2 

14af 


23.1a? 


23.61a? 


Be BsK^ 


2.853(2.816Jai 


2.1ai 


2.87a? 




BJ -» B^K 


0. 069(0. 061Jaf 


0.03ai 


0.13a? 




Be ^ -D tt" 


2.4(2. 0](xlO"*)a^ 
7.0(6. 0](xlO"'')ai 


- 


- 


- 




- 


- 


- 


BJ -> D+tj 


5.5(4.7](xl0 )a2 
3.1(2. 7](xl0"'')ai 


- 


- 






- 


- 


- 


B+ D+r]' 


1.7[1.5](xl0"'')a2 


- 


- 


- 




0.219(0. 185]ai 


- 


- 


- 


S+ — » D+D^ 


0. 261(0. 212]ai 
2.4(2.6](xl0"^)ai 


- 


- 




T-l-l- 7-1-1- 


- 


- 


- 


B+ ^ 


7.1[7.6](xl0"^)ai 


- 


- 




-De 


5.6(d.0](x10 ja2 








Bc+ ^ I^ + T) 


3.2[3.4](xl0-^)ai 








B+ ^ 


1.7[1.9](xl0-'5)ai 








B+ ^ D+D° 


0.0216(0.0227]ai 








B+ ^ D+L>*o 


0.0248[0.0250]ai 
0.779[0.648]ai 








II B+ ^ B+^» 


0.5ai 


0.54ai 


1.03ai 


5+ ^ B+p° 


0.967[0.752]ai 


0.7ai 


0.71ai 


1.28ai 


B+ B+Lj 


0.721[0.558]ai 








B+ B+T] 


3.99[3.30]ai 








B+ ^ B+v' 


0.054[0.045]ai 








B+ ^ B+K^ 


0.125[0.104]ai 








B+ ^ B+K" 


47.85[39.66]ai 


34ai 


35.3ai 




B+ ^ B+7^*o 


0.040[0.030]ai 










15.36[11.38]ai 


13ai 


13.1ai 









(O.Ollai + 0.011a2)^[(0.0097ai + 0.0097a2)^] 
(0.058ai + 0.064a2)^((0.052ai + 0.066a2)^] 
(0.428ai + 0.226a2H(0-482ai + 0.208a2f ] 
(2.27ai + 1.32a2H(2-47ai -f- 1.34a2)^] 




III 5+ 






B+ 




(0.438ai -H0.236a2)^ (0.47ai + 0.73a2)^ 


Bt 


^ VcDt 


(2.54ai -Hl.93a2)^ (2.59ai + 3.40a2f 



decays. It is thus expected that the dominant contri- 
bution to the Be total rate comes from the c induced 
decays. The more c induced Be VP and Be — > VV 
decay modes seem to deserve further consideration. 
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TABLE VI; Branching ratios (in %) of the exclusive non-leptonic Be decays at the fixed choice of Wilson coefficients: af (aj) = 
1.20 (1.14) and a2 ('^2) = —0.317 (—0.20) relevant for the nonleptonic decays of the c (b) quark. For the lifetime of the Be we 
take t{Bc) = 0.46 ps. 



Class Mode 


Lin [HO] 


[8,9] 






[11 




[5] 




[13] 


[40] 




[10] 




B+ 


D%+ 


4.3[4.3](xl0-'') 


- 


- 


- 


- 




- 




- 


- 




- 




B+ 




1.3[l.l](xl0-*) 


- 


- 


- 


- 




- 




- 


- 




- 




B+ 




3.5[3.1](xl0-'') 


- 


- 


- 


- 




- 




- 


- 




- 




B+ 


^ D°K*+ 


6.8[5.8](xl0"'') 


- 


- 


- 


- 




- 




- 


- 




- 




B+ 




0.091(0.116] 


0.085 


0.094 


0.13 


0.19 




0.20 




0.025 


- 




0.18 




Bt 




0.257(0.324] 


0.21 


0.24 


0.30 


0.45 




0.42 




0.067 


- 




0.49 




Bt 




0.0074(0.0094] 


0.0075 


0.0075 


0.013 


0.015 




0.013 




0.002 


- 




0.014 


I 


Bt 




0.013(0.017] 


0.011 


0.013 


0.021 


0.025 




0.020 




0.004 


- 




0.025 




B+ 




0.157(0.131] 


0.10 


0.11 


0.15 


0.20 




1.06 




0.19 


0.373 




0.32 




Bt 


^bV 


0.195(0.152] 


0.13 


0.14 


0.19 


0.20 




0.96 




0.15 


0.527 




0.59 




Bt 


-*B''K+ 


0.013(0.011] 


0.009 


0.010 


- 


0.015 




0.07 




0.014 


0.027 




0.025 




Bt 




0.0042(0.0032] 


0.004 


0.0039 


- 


0.0048 




0.015 




0.003 


0.023 




0.018 




Bt 




3.723(3.697] 


2.52 


3.51 


3.42 


3.9 




16.4 




3.01 


5.309 




5.75 




Bt 


- B°P+ 


2.561(2.338] 


1.41 


2.34 


2.33 


2.3 




7.2 




1.34 


6.265 




4.41 




Bt 


- B°K+ 


0.287(0.284] 


0.21 


0.29 


- 


0.29 




1.06 




0.21 


0.367 




0.41 




Bt 


^ B°K'+ 


0.0069(0.0061] 


0.003 


0.013 


- 


0.011 




- 




0.0043 


0.165 




- 




Bt 




6.7(5.6](xl0~'') 


- 


- 


- 


- 




- 




- 


- 




- 




Bt 


~*D+p° 


2.0[1.7](xl0"'') 


- 


- 


- 


- 




- 




- 


- 




- 




Bt 


D + LJ 


1.5[1.3](xl0"'') 


- 


- 


- 


- 




- 




- 


- 




- 




Bt 


D+T] 


8.7[7.6](xl0"'^) 


- 


- 


- 


- 




- 




- 


- 




- 




Bt 




4.8[4.2](xl0"*) 


- 


- 


- 


- 




- 




- 


- 




- 




Bt 


-*D+D° 


6.1[5.2](xl0-'') 


- 


- 


- 


3.3 X 10" 


-3 


5.3 X 10" 


-3 


4.1 X 10"* 


- 




1.8 X 10"=* 




Bt 


^ D+D*° 


7.3[5.9](xl0-'') 


- 


- 


- 


3.8 x 10" 


-3 


7.5 X 10" 


-3 


3.6 X 10""* 


- 




1.9 X 10"^ 




Bt 


- Dt^° 


6.7[7.3](xl0-*) 


- 


- 


- 


- 




- 




- 


- 




- 




Bt 


- Dtp° 


2.0[2.1](xl0"'^) 


- 


- 


- 


- 




- 




- 


- 




- 




Bt 


^Dtuj 


1.6(1. 7](xl0-'^) 


- 


- 


- 


- 




- 




- 


- 




- 




Bt 


^ Dtv 


9.0[9.5](xl0-*) 


- 


- 


- 


- 




- 




- 


- 




- 




Bt 


- Dtv' 


4.8[5.3](xl0-*) 


- 


- 


- 


- 




- 




- 


- 




- 




Bt 


- DtD'> 


6.0[6.3](xl0~^) 


- 


- 


- 


2.1 X 10" 


-4 


4.8 X 10" 


-4 


2.7 X 10"^ 


- 




9.3 X lO"'^ 




Bt 


-> DtD*^ 


6.9[7.0](xl0"^) 


- 


- 


- 


2.4 X 10" 


-4 


7.1 X 10" 


-4 


2.5 X lO"'^ 


- 




9.7 X lO"'^ 


II 


Bt 


^B + 7V° 


0.0055(0.0046] 


0.004 


0.0038 


0.007 


0.007 




0.037 




- 


4.6 X 10" 


-5 


0.011 




Bt 


-^B+p° 


0.0068(0.0053] 


0.005 


0.0050 


0.009 


0.0071 




0.034 




- 


6.5 X 10" 


-5 


0.020 




Bt 


B + LU 


0.0051(0.0039] 


- 


- 


- 


- 




- 




- 


5.8 X 10" 


-5 


- 




Bt 


B+T] 


0.028(0.023] 
o.o[o.z\[X W ) 


- 


- 


- 


- 




- 




- 


1.6 X 10" 


-4 


- 




Bt 


^B+v' 




















-6 






Bt 


^B+K° 


8.8[7.3](xl0-'') 


















6.5 X 10" 


-6 






Bt 


^B+K" 


0.336(0.279] 
2.8(2.1](xl0-'') 


0.24 


0.25 




0.38 




1.98 






0.0022 




0.66 




Bt 


^ B+K*" 


















5.5 X 10" 


'6 






Bt 


-> B+K*" 


0.108(0.080] 


0.09 


0.093 




0.11 




0.43 






0.0018 




0.47 




Bt 


^ D+D' 


7.5(5.8](xl0"'') 








3.1 X 10" 


-5 


3.2 X 10" 












III 


Bt 


^ DtD° 


2.0[1.5](xl0"'') 








7.4 X 10" 


-4 


6.6 X 10" 


-4 












Bt 




0.014(0.018] 




0.014 


0.010 


0.019 




0.032 




0.0055 






0.0012 




Bt 


^ VcDt 


0.378(0.454] 




0.44 


0.35 


0.44 




0.86 




0.51 






0.056 
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